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Polymers in a melt may be subject to topological constraints, as in the example of unlinked polymer 
rings. How to do statistical mechanics in the presence of such constraints remains a fundamental 
open problem. We study the effect of topological constraints on a melt of directed polymers, using 
simulations of a simple quasi-2D model. We find that fixing the global topology of the melt to 
be trivial changes the polymer conformations drastically. Polymers of length L wander in the 
transverse direction only by a distance of order (InL)*" with ~ 1.5. This is strongly suppressed in 
comparison with the Brownian scaling which holds in the absence of the topological constraint. 
It is also much smaller than the predictions of standard heuristic approaches — in particular the 
lG'I of a mean-held-like ‘array of obstacles’ model — so our results present a sharp challenge to 
theory. Dynamics are also strongly affected by the constraints, and a tagged monomer in an inhnite 
system performs logarithmically slow subdiffusion in the transverse direction. To cast light on the 
suppression of the strands’ wandering, we analyse the topological complexity of subregions of the 
melt: the complexity is also logarithmically small, and is related to the wandering by a power law. 
We comment on insights the results give for 3D melts, directed and non-directed. 


The fact that polymer chains cannot pass through each 
other is the crucial factor in their dynamics, underlying 
for example the reptation picture m, and in various 
situations also determines their equilibrium state. Two 
salient examples are a single ring polymer and a melt 
of rings that do not knot or link. In such cases the no¬ 
crossing condition sets topological constraints that are 
inherently non-local. The statistical mechanics of such 
systems is a tremendous theoretical challenge, for which 
no systematic theoretical tools are presently available. 
Heuristic approaches [IHID] and ever-growing numeri¬ 
cal simulations mHiT] have provided substantial insight, 
but even basic issues — such as the size of a single ring 
polymer in a melt or the degree to which different rings 
mix — are not resolved. Ring polymer melts have re¬ 
ceived considerable attention as models of chromosome 
arrangement in the nucleus CHI: and experiments on ring 
melts have revealed unique rheological properties m- 
Additionally, dense systems of open chains may be sub¬ 
ject to effective topological constraints on intermediate 
timescales, yielding very slow relaxation and long-lived 
‘pseudoequilibrium’ states with less entanglement than 
at equilibrium [2011^ . 


state in which all polymers are straight lines. The melt 
is endowed with Monte-Carlo dynamics that respect this 
constraint (i.e., respect the fact that the polymers cannot 
pass through each other). Mathematically, the polymers 
form a trivial ‘braid’. The statistical properties of ran¬ 
dom braids have been studied extensively in order to shed 
light on polymer topology [22H26] , but the dynamical and 
conformational properties of a topologically constrained 
melt have not been investigated. 

Our model is extremely tractable computationally, so 
we are able to obtain precise results for its universal prop¬ 
erties: these turn out to be surprising in the light of cur¬ 
rent theoretical ideas. The model is also simple enough 
to allow hope of analytical progress. 

A key feature of the model is that it allows comparison 
with the predictions of standard theoretical approaches, 
shedding light on the validity of ideas that are more gen¬ 
eral than the directed case: for example the idea of mod¬ 
elling topological effects in a melt using toy models of a 
single polymer in an array of obstacles 13 El HQ], or the 
use of Flory-like arguments [3i]. These approaches are 
widespread but hard to justify a priori. It is important 
to find ways to confront them with precise results from 


The aim of this paper is to study the simplest possi¬ 
ble (but genuinely many-body) model for a topologically 
constrained melt. Physically, this model describes di¬ 
rected polymers in quasi-2D, i.e. in a slab geometry, 
but with the positions of the polymers projected onto 
the plane to give a 2D lattice model. The remnants of 
three-dimensionality are the fact that the polymers can 
cross and the crucial distinction between over and un¬ 
dercrossings (Fig. [^. The endpoints of the polymers are 
held fixed, and the entire melt is constrained to be topo¬ 
logically trivial: that is, continuously deformable to the 



FIG. 1. (a) Topologically trivial configuration in a small sys¬ 
tem. (b) Monte Carlo move types. There are six variants of 
the move on the lower right, and two of that on upper right. 
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a genuinely interacting many-polymer model. 

The present model may also capture the universal be¬ 
haviour of some realistic situations. The most striking 
features of our results are expected to extend to systems 
of directed polymers in 3D, so the results for dynamical 
behaviour may be relevant to relaxation and equilibra¬ 
tion in polymer brushes [23128] . 

In both the present model and 3D ring melts, topolog¬ 
ical constraints reduce the extension of individual poly¬ 
mers: entropy dictates that polymers ‘hide’ from each 
other, as more extended configurations are more likely to 
be entangled. We will soon see that for directed polymers 
this effect is almost as strong as it could possibly be. 

Model. Take a number of strands, each of height 
Ly and directed in the y direction, see Fig. [^a). At in¬ 
teger j/-values the polymers lie at integer cc-values, and 
all lattice points are occupied. Between y and y -I- 1, a 
given polymer may be vertical, or two neighbours may 
cross. We take periodic boundary conditions in the x di¬ 
rection, and fix the positions of the endpoints at y = 0 
and y = Ly so the two ends of a given polymer have the 
same a:-coordinate. Finally we enforce topological trivi¬ 
ality: the allowed configurations are those which can be 
deformed to the configuration with straight vertical poly¬ 
mers. Mathematically, each configuration C defines an 
element g{C) of the braid group [22], and allowed con¬ 
figurations are those in which this is the trivial element 
‘1’. The partition function Z is the equally weighted sum 
over allowed configurations {6 is the Kronecker delta): 

^9{C},1 ( 1 ) 

configs C 

In practice, to fully define the model we need a Monte 
Carlo scheme which samples only the topologically triv¬ 
ial configurations. For this we use the moves shown in 
Fig.0b) , which have a simple relationship with the defin¬ 
ing relations of the braid group |30j . and which form a 
complete set of moves. Intuitively, any local rearrange¬ 
ment can be decomposed into: creation/annihilation of 
pairs of crossings; motion of crossings; and motion of one 



FIG. 2. A subregion of a topologically trivial melt (taken 
from a system of size = 512, Ly = 128). 
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FIG. 3. Main panel: transverse wandering D of a strand at 
its midpoint y = Lyl2, as a function of the length Ly of the 
strands (note lin-log scale). Red line: fit to the logarithmic 
form in Eq. [^ Inset: effective exponent = dlnD/dlnLy, 
see text, with fit derived from Eq. [^ 


strand over/under a crossing between two others. These 
are precisely the moves in Fig. 0b)- Fig.l gives an idea 
of what a subregion of the melt looks like. Note that 
the model has a fixed monomer density: this eliminates 
finite-size effects due to fluctuations of the density mode, 
which is irrelevant to long-distance behaviour. 

We work in the limit where L^, the number of strands, 
is much greater than the typical wandering of the strands. 
In fact we enforce a stricter criterion: we ensure L^ is 
large enough that the results are essentially those of the 
La; —> oo limit. We find that this can be achieved using 
modest L^, which is unsurprising given that the wander¬ 
ing is much smaller than Ly. Below we take Ly rang¬ 
ing up to Ly = 1200, and L^ ranging up to L^ = 100 
(larger for some Ly). Appendix [A] gives further details of 
simulations, including basic checks of equilibration and 
convergence in L^. 

Results. The wandering Di{y) is defined as the trans¬ 
verse displacement of the fth strand at height y. We 
denote the root mean square transverse displacement by 
D{y). In the absence of the topological constraint in 
Eq. [^ each strand essentially performs a random walk 
constrained to return to Xy^o when y = Ly\ thus it is 
clear (and we have checked, App. [0 that in the uncon¬ 
strained case D{Ly/2) ~ a/L^. 

As noted above, the topological constraint will reduce 
the polymers’ wandering. Fig. [^ quantifies this using 
the r.m.s. wandering at the midpoint of the strands, 
D{Ly/2), plotted against Ly. The fit is 


c = 1.49(3), (2) 
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FIG. 4. Main panel: r.m.s. wandering D{y) versus position 
y along the strand, for various values of the strand length 
Ly. Blue line: fit to H = A' {Iny /Iq)‘’ for 10 < t/ < 100, 
giving C' = 1.54(15), A' = 0.204(10), I'o = 0.146(10). Lower 
inset: same data on log-log plot. Upper inset: probability 
distribution for wandering of a strand, P{Di), rescaled by 
standard deviation D. 

with A = 0.26(2), lo = 0.91(9). 

This logarithmic form is unexpected. However, a more 
conventional power law fit D oc Ly leads to much worse 
results, or to an exponent equal to zero within error bars 
if subleading corrections are included [32]. The inset to 
Fig. [3] shows an effective finite-size wandering exponent 
defined by VeS = dlnU/dlniy. This drifts downwards, 
as expected for the logarithmic form, according to which 
r'eff —>■ 0 as Ly —oo. The ‘array-of-obstacles’ prediction 
discussed below, z/ = 1/4, is clearly ruled out. See App.|A| 
for further discussion of fits. 

Further evidence for the logarithmic behaviour comes 
from the wandering D{jj) in the distinct regime y Ly. 
Since no scaling theory exists for this problem, it is 
not guaranteed a priori that the behaviour for y Ly 
and y ~ Fy/2 will be similar, but this turns out to 
be the case. See Fig. [^ where the data fits well to 
D{y) = (iny/Vof with C' = 1.54(15). Note the strik¬ 
ing agreement between the independently-determined ex¬ 
ponents and C'. The lower inset to Fig. shows the 
same data on a log-log scale: it is clear from the curva¬ 
ture that a power law would fit only over a very narrow 
range of scales. The result for D{y) highlights the fact 
that the properties of a finite-sized subsystem (of height 
y) are strongly affected by the global topological con¬ 
straint even in the limit Ly —> oo. 

In addition to the r.m.s. displacement of a strand we 
may consider the full probability distribution. Fig. [^ 
(upper inset) shows this for the displacement at y = 
Ly/2. The data collapses beautifully after rescaling by 
the r.m.s. value. The distribution is not quite Gaus¬ 


sian (App. 0. 

Correlations. Correlations between the displacements 
of different strands decay exponentially when their sep¬ 
aration is much larger than the wandering. Specifically, 
let Cd{x) = {DiDi+x), where Di is the transverse dis¬ 
placement of the zth strand at its midpoint. At large 
X, Cd{x) ^ with a correlation length ^(Ly) that 

grows in a roughly similar manner to D (App. 0. The 
two-point function Cx(x) for the density of crossings (for 
plaquettes at y = Ly/2 separated by a distance x) decays 
exponentially with period 2 oscillations and a correlation 
length of less than two lattice spacings. 

Dynamics and logarithmic subdiffusion. The timescale 
T for relaxation of the melt is extracted from the Monte 
Carlo time series. It is independent of L^ for large L^, 
but depends nontrivially on the length Ly of the strands: 

r ~ L;, z = 2.60(6). (3) 

This exponent describes the equilibration of the en¬ 
tire system within the topologically constrained space of 
states. The transverse motion of a tagged monomer in 
an infinite system is more interesting. By Eq. we ex¬ 
pect that motion of the monomer by a distance x involves 
rearrangements of segments of height y ~ exp(x/A)^A^ 
and a time t which is a power law in y. This implies that 
the tagged monomer subdiffuses logarithmically slowly: 

(x^) - {lntf‘^ ■ (4) 

Similarly, the decorrelation time of a subregion of the 
melt of size x x y is exponentially large in x. 

Eq. 1^ describes a monomer inside the topologically 
trivial melt. It is interesting to ask about the dynam¬ 
ics of a monomer for other choices of the initial state; 
for example an equilibrated state of the topologically 
unconstrained problem. A naive guess might be that the 
increased local entanglement in such a state will slow the 
motion even further, but this needs investigation. These 
issues are related to relaxation in polymer brushes m, 
in which the polymers are tethered at one end and are 
directed (on large scales) for high surface fraction [2711^ . 

Topological complexity. The present model allows for 
a clean definition of the topological complexity of a sub- 
region (the full system is of course topologically trivial). 
Examining this complexity illuminates the drastic sup¬ 
pression of the wandering. We take the subregion to 
be the bottom half of the braid, y < Ly/2. Eixing the 
strands’ endpoints at y = 0 and y = Ly/2 gives the half¬ 
braid a well-defined topology which cannot be changed 
by allowed moves in the interior. Allowed moves can how¬ 
ever reduce the total number of crossings. Let A^min be 
the minimal value to which we can reduce this number, 
and define the complexity per strand as 

C = 2A^iiii „/Lx 


(5) 
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FIG. 5. Main panel: complexity (per strand) for the half¬ 
braid as a function of Ly, fitted to the form in Eq. with 
77 = 2.6(2). Inset: Wandering D{Ly/2) plotted versus the 
complexity C{Lyl2) and fitted to a power law, Eq. 

C is the average number of crossings encountered (i.e. 
steps to the right or left taken) by a single strand in 
the reduced configuration. C is finite as ^ oo. Note 
that reducing the half-braid does not change Di{Ly/2): 
strands in the reduced half-braid take many fewer steps, 
but wander by the same total distance. 

We compute C by simulated annealing. Starting with 
an equilibrated braid, we extract the lower half and sub¬ 
ject it to a modified Monte Carlo dynamics with an en¬ 
ergy penalty for crossings. The temperature is grad¬ 
ually lowered until the system finds its ‘ground state’. 
We do not encounter problems finding the ground state 
(App.g, perhaps because C is modest. 

The main panel of Fig. shows C plotted against Ly, 
and the inset shows the r.m.s. wandering D(Ly/2) plot¬ 
ted against C. Strikingly, the wandering has a clean 
power law dependence on the topological complexity: 

D(xC‘^, 0 = 0.618(2). ( 6 ) 

This implies, for consistency with Eq. 

C(x(\n^\ , ?7 =-~ 2.41(6). (7) 

\ lo J a 

This is indeed compatible with the results in Fig. For 
comparison, a braid configuration of height Tj ,/2 chosen 
uniformly from the set of all such braids has a typical 
complexity of order Ly |23L I24j . For a fixed finite number 
of strands, a sub-braid of sufficiently large height y is 
believed to have a complexity of order ^ (THl HH US] , 
but this is a different (‘quasi-lD’) limit (App. 0 . 

The suppression of wandering may therefore be viewed 
as a consequence of the drastic suppression of complexity. 
Interestingly, the geometry of the strands in the reduced 
half-braid is more conventional than in the unreduced 


half-braid: the wandering I? of a strand has a power law 
dependence on the average number of steps, C. The ex¬ 
ponent a in Eq. I^is greater than 1/2, indicating positive 
correlations between steps in the reduced half-braid. By 
contrast, the number of steps per strand in the unreduced 
half-braid is much larger, 0{Ly), and there are strong 
negative correlations between steps. 

Comparison with standard ideas. Many approaches to 
topologically constrained melts rely either on simplifying 
the problem to a single strand in an array of obstacles 
which represent the other polymers nnnniEa , or on 
Flory-like free energy arguments [THg. These ideas have 
for example been used to argue that rings in a 3D melt 
fold up into compact tree-like structures [3]. However 
both approaches are uncontrolled approximations which 
must be tested against data. Eor 3D ring melts this is 
challenging because of large finite-size effects [HI [33] . 

Here, we can make a quantitative comparison with the 
natural array of obstacles model for the directed case, 
which describes a single fluctuating strand in an array of 
straight vertical strands (App.[D|). This predicts D ~ L)) 

with V = 1/4, and C ~ |S|, contrary to our results. 

The wandering distribution also differs (App. 0. The 
exponent 2 ; on the other hand is roughly compatible with 
the z = 5/2 of the array of obstacles model [101 (though, 
by the reasoning preceding Eq. [^ the transverse diffusion 
in the topologically constrained ensemble will be much 
faster for the array of obstacles). Our results show that, 
for directed polymers, the behaviour of the true melt is 
very different from the array of obstacles model. 

Flory-like estimates are sensitive to the assumed form 
of the entropic cost of not being entangled [3f], which 
is hard to control. Here one can easily show that 2- 
strand interactions alone are not enough to explain v = 0. 
See App. [^ for further discussion. The approach does 
however support the expectation that wandering is at 
least as strongly suppressed in 2-l-lD as in 1-l-lD. 

Throwing out 3-strand moves. Though much simpler 
than a melt of 3D rings, the present model is still a 
formidable challenge analytically. One may also con¬ 
sider a reduced model based on the ‘locally free group’ 
[23l[24|, a simplification of the braid group. This means 
imposing the stronger constraint that the melt be de¬ 
formable to the straight-line state without using the 3- 
strand moves of Fig. [^ This is a drastic simplifica¬ 
tion, and no longer faithful to the topology of directed 
melts. Nevertheless preliminary simulations suggest that 
behaviour for D remains qualitatively similar, with a re¬ 
duced ^ ^ 1 (App. H- 

Future directions. We believe that the crucial features 
of the present model, including the fact that the wan¬ 
dering is logarithmic (though not necessarily the value 
of C) will carry over to the 3D directed case. This is be¬ 
cause the number of other strands encountered by a given 
strand grows faster with D in 3D than in 2D, indicating a 
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stronger entropic penalty for wandering. This conjecture 
must be examined numerically. Another natural next 
step is to investigate the dynamics of the present model 
when the endpoints of the chains free to move, so that 
the topology of the melt can slowly relax. It would be 
interesting to know whether the transverse motion of the 
monomers remains logarithmically slow even in the final 
equilibrated state. If we start from an unentangled con¬ 
figuration, even the static properties may remain similar 
to those discussed here for a very long time. 

We have seen that for a topologically constrained en¬ 
semble of directed polymers, the exponent governing the 
chains’ extension takes its minimal possible value, z/ = 0, 
with logarithmic corrections. One might wonder whether 
in a 3D ring melt the radius of gyration is also governed 
by logarithmic corrections to the minimal exponent value 
(i/ = 1/3). Conceivably, such logarithms might partially 
explain the slow saturation observed for v in this case. 

A fundamental question is whether there exists a real- 
space renormalisation group treatment for topologically 
constrained polymers (App. 0. The present model may 
be simple enough to offer hope of this. 
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Appendix A: Further details of simulations 

The appendices include: further details of simulations, 
fits etc. (this section); a discussion of the effect of 3- 
strand moves (App. 0; comparison with the topolog¬ 
ically unconstrained case (App. 0 ; and further infor¬ 
mation about the heuristic approaches mentioned in the 
main text (App. 0. 

Equilibration. An important step is to ensure the equi¬ 
libration of the sample. To begin with, we calculate a re¬ 
laxation time Ly) from the time series of the r.m.s. 

wandering I? in a standard way, by fitting the autocor¬ 
relation function to an exponential. For the Ly values 
considered, T{Lx,Ly) rapidly approaches a constant in¬ 
dependent of Lx (on a scale much less than 100). The 
Ly-dependence of the autocorrelation time is plotted in 
the inset to Fig. 

For a stringent check on the equilibration of our sam¬ 
ples, at each system size we run simulations starting from 
two very different (topologically trivial) initial states. 



FIG. 6. Main panel: correlation lengths ^ for (1) correlation 
Cd between the wandering of two strands, and (2) correlator 
Cx of the density of crossings. The fit (red line) for (1) is to 
= a(lnLj,/Zo)*' with a = 0.14(8), lo = 2.1(11), b = 1.8(2). 
Inset: correlation time as a function of Ly, fitted to t oc Ly 
with 2 = 2.60(6). 


The first is the conhguration where all polymers are 
straight vertical lines. The second is constructed by gen¬ 
erating a lower half with crossings at random and taking 
the upper half to be the mirror image of the lower one. 
This gives a configuration which is globally trivial but 
locally highly entangled. The two cases represent oppo¬ 
site extremes both for the wandering at y = Ly/2 (which 
is exactly zero in the first case and in the sec¬ 

ond) and for the half-braid complexity C (zero and 0{Ly) 
respectively). Nevertheless they converge to the same 
equilibrated value of the wandering under the Monte 
Garlo dynamics. After 20 autocorrelation times, results 
from straight lines and maximally disordered samples are 
identical to within error bars: the inset to Fig. shows 
(Da — D)/D, where Da is the wandering for samples with 
different initial conhgurations a = 1,2 (straight lines or 
disordered) and D is the average of both. We therefore 
start collecting data after 20t Monte Garlo sweeps. 

Convergence in Lx. We examined the dependence of 
D and other quantities on the horizontal system size Lx 
in order to ensure that the results are those of the large 
Lx limit. We studied the nature of this convergence ex¬ 
tensively for Ly = 6, 32 and 100. We see exponential 
convergence in Lx with a very small characteristic length 
and signs of oscillations within the envelope. The main 
panel of Fig.j^shows the wandering for Ly = 100 plotted 
as a function of Lx', the red line is an exponential fit with 
a characteristic length Zq = 2.8 ± 1.1. For larger Ly val¬ 
ues, we checked that varying Lx did not change the value 
of D, within error bars. (For most samples we considered 
Lx < 100, and for Ly = 800 we considered Lx < 200.) 
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Correlation functions also yield short correlation 
lengths, as discussed in the text and shown in Fig.[^ 

Non-Gaussianity of P{Di). The probability dis¬ 
tribution P{Di) in Fig. 4 of the main text is close 
to, but measurably different from, a Gaussian. To 
see this we examine the difference of moment ra¬ 
tios M.k = {\Di\^) / from the Gaussian value: 

A4fc = M-k — ^Gaussian^ pg^. comparison, we have also 
obtained these moments for the toy model of a single 
strand in an array of obstacles (see below), by a separate 
simulation. 

For A^ 3 / 2 , an extrapolation to Ly = oo gives —0.021(1) 
for the full model and —0.0363(3) for the array of ob¬ 
stacles, indicating that the universal scaling function is 
different in the two cases. The error bars are statistical 
errors in extrapolations of the form M. = A + BL~^; sys¬ 
tematic errors may be larger. For the fourth cumulant 
the extrapolation to Ly = oo is more difficult but we find 
M 4 ~ 0.2 for the full model and M 4 ~ 0.5 for the toy 
model. 

Computation of topological complexity of half-braid. To 
check the simulated annealing protocol described in the 
text, we reduce the same half-braid multiple times and 
confirm that after each iteration the same number of 
crossings remains. We also check that the protocol suc¬ 
ceeds in reducing the full braid to the straight-line con¬ 
figuration. Finally we check that reducing the rate at 
which the temperature is decreased does not change the 
results. For illustrative purposes. Fig. shows a small 
system before and after reduction. 

It has been conjectured that for a fixed finite number 
of strands, in the limit of large Ly, subregions of height 
?/ of a trivial braid have a complexity proportional to 



FIG. 7. Main Panel: Wandering D as a function of for 
samples of size x 100. Inset: Relative error in wandering 
for the two different initial configurations, {Da — D)ID, as a 
function oi Ly (a = 1, 2). 



FIG. 8. The reduction of a half-braid. (For illustrative pur¬ 
poses only: results in text use larger values.) 


for sufficiently large y [m Ea [23]. Our result that C is 
only logarithmically large is surprising in the light of this, 
but the results are not in contradiction: the former is for 
the ‘quasi-lD’ situation, in which Lx is fixed as Ly —> 00 , 
whereas we consider the 2D situation in which Lx is much 
larger than the scale of the wandering of the strands. It 
would be interesting to understand the crossover between 
the two limits. 

Further checks on logarithmic fits. In the main text we 
fitted the wandering (Figs. 3, 4 of the main text) to the 
forms D{Ly/2) = A{\iiLy/lo)‘' and D{y) = A'{\ny/I q)'^ , 
finding remarkable agreement between the exponents C 
and Cf. Here we check the extent to which finite-size 
corrections to these forms could affect the values of C, and 
C,'. (We have already noted in the text that power-law 
fits are much less convincing than logarithmic ones.) 

We have tried various possible forms for sub¬ 
leading corrections. Fits with C = 1, such as 

Hln(Ly//o) + H/(lnLj,//o) or A\D.{Ly/lQ) + B/Ly, are 
poor. For a more stringent test we consider 


and the analogue for D{y). The resulting exponent values 
are C = 1.49(7) and C,' = 1.59(22), where the error bars 
are calculated using the bootstrap method. For the 
range used for the fit is Ly > 40 (we have checked that 
the exponent is stable when this value is increased) and 
for C', where we must ensure y <C Ly, it is 10 < y < 100. 

The fact that the values of C and Q' are stable under the 
addition of subleading terms to the fit, and in particular 
the agreement between C, and C,', gives us confidence that 
these exponents are indeed the same, and that they are 
close to 1.5 rather than being equal to any integer value. 
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FIG. 9. Comparison of the wandering D{Ly/2) for the full 
model (same data as in the main text, black) and the model 
without the three-strand moves (red). Fit described in text. 


Appendix B: Simulations without 3-strand moves 


FIG. 10. Main Panel: Wandering as a function of system size 
in the topologically unconstrained problem, fitted to a power- 
law with exponent 0.514(2). Inset: Wandering as a function 
of y for several system sizes: straight line is a power law with 
Brownian exponent 0.5. 


We have performed limited simulations without using the 
3-strand moves of Fig. 1 of the main text, starting with 
the vertical line initial condition. This corresponds to 
a partition function in which the polymer configurations 
satisfy a stronger constraint. This constraint does not 
have as natural an interpretation in terms of topology of 
strands: it means that strands are forbidden from mov¬ 
ing over crossings between other pairs. However it has 
been suggested as a natural simplification of the alge¬ 
braic structure of the problem, equivalent to replacing 
the braid group with the locally free group [231123] • 

Fig. [^compares D{Ly/2) for the full model discussed in 
the main text with that for the model without 3-strand 
moves. A fit to the form A”{lnLy/lQ)‘^ for the latter 
gives A" = 0.5(3), = 0.9(1.1), and C" = 1-0(3). 


Appendix C: Unconstrained model 


analytically, giving a useful check on our algorithms. 

We also checked the case where the topology-changing 
moves are allowed but the endpoints are kept fixed. The 
corresponding partition function is a sum over config¬ 
urations of any topology but with fixed endpoints. The 
wandering remains Brownian and the correlators Cjy and 
Cx still decay exponentially. The wandering D{Ly/2) 
for this topologically unconstrained case is plotted in the 
main panel of Fig. 10 as a function of Ly. The results 
shown are for samples with = 200 (convergence to 
the limit La, —> oo is again fast). As expected, the wan¬ 
dering is compatible with Brownian behaviour D ^ Ly 
with 1 / = 1/2 (a fit for Ly > 16 gives exponent z/ ~ 0.51; 
the exponent decreases towards 0.5 when dropping small 
system sizes). As a further check the inset shows be¬ 
haviour for y Ly, again giving results compatible with 
Brownian behaviour. 


As basic check of our algorithms, and to confirm that 
the unusual features found in the text (e.g. the dras¬ 
tic suppression of wandering) are due to the topological 
constraint, we have characterised the model without this 
constraint. As expected it shows simple Brownian wan¬ 
dering as a function of the height y (not to be confused 
with dynamics in time t), and exponentially decaying 
correlation functions. 

To do this we allow an additional move that exchanges 
undercrossings and overcrossings. If we also relax the 
boundary conditions at the top of the sample (so that 
a strand is allowed to have different a:-coordinates at 
y = 0 and y = Ly), the model becomes easily solvable 


Appendix D: Heuristic approaches 


Array of obstacles. Fig. 11 illustrates this toy model. 
There are fixed vertical strands at integer values of the 
x-coordinate, and the mobile strand winds in and out. 
The endpoints of the mobile strand at y = 0, y = Ly are 
fixed at x = 1/2, and the global topology is constrained 
to be trivial (the mobile strand can be deformed to ver¬ 
tical without passing through any background strands). 
Variants of this model have been considered many times 
in the literature IMZllM]. Remarkably, the mathemat¬ 
ical structure of this toy model is identical to that of 
the array of obstacles model for the 3D undirected case. 













FIG. 11. Toy ‘array of obstacles’ model. 

despite the different physical interpretations. The exact 
exponents quoted in the text {v = 1/4 and C ^ 
come from an elegant mapping to random walks on the 
Cayley tree [S]. 

Entropic arguments. A naive approach to our model is 
to minimise a free energy per strand, F = T/onfinement + 
Ttopoiogicai, as a function of the typical scale of wandering, 
D. The first term is the cost of confining a free directed 
strand in a Z?~sized ‘box’: ^confinement Ly/D^. The 
second is the free energy cost of not being entangled — 
i.e. of of obeying the constraint that the melt is topolog¬ 
ically trivial — and is harder to estimate. The constraint 
of topological triviality can be regarded as a union of an 
infinite number of ‘k-strand’ constraints, where by a k- 
strand constraint we mean the requirement that a given 
group of k strands is not entangled. 

First let us confirm that the two-strand constraints 
alone are far too weak to generate logarithmic confine¬ 
ment (two-strand constraints should not be confused 
with the two-strand moves discussed in Sec. II above). 
If d is the number of transverse dimensions (d = 1 in 
the model simulated), our chosen strand has the oppor¬ 
tunity to wind around 0{D‘^) other strands, each one 
0{Ly/D‘^) times. Viewing y as time, the winding num¬ 
ber of a given pair performs a ID random walk with 
0{Ly/D‘^) steps, which must return to the origin. This 
gives F 2 _strand ~ ln{Ly /D‘^). Minimising F gives 

D ^ L'^ (neglecting logarithms) with v = 1/(2-|-d). This 
confirms that the two-strand constraints are too weak to 
reproduce what we see: genuinely multi-strand entangle¬ 
ment is crucial. The above result does however support 
the expectation that wandering is at least as strongly 
suppressed in 3D as in 2D, as a result of the larger num¬ 
ber of other strands encountered for a given typical scale 
of wandering. 

Attempting to go beyond the two strand constraints il¬ 
lustrates the fact that the results will be sensitive to the 
approximate form assumed for the free energy. We may 


imagine successively imposing the 2-strand constraints, 
then the 3-strand constraints that are independent of 
them, then the additional 4-strand constraints etc. At 
each stage we pay a free energy Ffe.gtrand for forbidding 
configurations with ‘fc-partite’ entanglement. The crud¬ 
est approximation is to approximate Ffe.gtrand as a sum 
of independent terms, one for each group of fc-strands 
which includes the chosen strand and which are not spa¬ 
tially disjoint. (Note that two distinct /c-strand groups 
may share some smaller subset of their strands.) This 
gives Ffc.strand ^ fk (for large D and d = 1), where 
fk is the free energy cost of forbidding fc-partite entan¬ 
glement for a single fc-strand group. It is tempting to 
imagine that a conservative approximation is to take fk 
to be at least of order one. Then, if we truncate the free 
energy at some fixed k, we obtain an exponent Vk which 
indeed tends to zero as fc —)■ oo. 

This suggests that interactions between arbitrarily 
large numbers of strands are necessary to explain the ob¬ 
served fact that iz = 0. However, the argument may be 
misleading; treating the groups as independent is an un¬ 
controlled approximation. The array of obstacles model 
shows the danger. Naively we might attempt a similar 
argument there (fc-strand entanglement can arise even if 
only one strand is mobile) but we know that in that case n 
is not equal to zero. It is an interesting question whether 
in the full model fc-strand constraints for arbitrarily large 
k are necessary to obtain v = 0. 

The renormalisation group. One of the most funda¬ 
mental questions is whether there exists a real-space 
renormalisation group approach to topologically con¬ 
strained polymer ensembles (see e.g. Ref. |35] for the 
case without constraints). Our result that D is only log¬ 
arithmically large in Ly suggests that the appropriate 
protocol here would be to rescale the vertical coordinate 
but not the horizontal one, i.e. not to decimate strands. 
One must of course find a way to deal with the nonlo¬ 
cality of the topological constraints. One possibility is 
not to directly coarse-grain the 2D configuration, but in¬ 
stead to focus on the configurations in 2-l-lD spacetime, 
with the same time evolution as the Monte Carlo simula¬ 
tions. Spacetime configurations may be viewed in terms 
of ‘worldlines’ of crossings. The potential simplification is 
that the global topological constraints are then encoded 
in local constraints on the worldline configurations. If 
the problem is simplified by dropping the three strand 
moves from the dynamics (as discussed in the text) we 
obtain a 2-l-lD statistical mechanics problem with a con¬ 
cise definition. This is a multi-layer loop model, with 
oriented loops, in which loops from adjacent layers can¬ 
not cross. It is conceivable that this model is tractable 
by field theory techniques. 
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